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The dynamics of a binary system with two spinning components on an eccentric orbit is studied, 
with the inclusion of the spin-spin interaction terms appearing at the second post-Newtonian order. 
A generalized true anomaly parametrization properly describes the radial component of the motion. 
The average over one radial period of the magnitude of the orbital angular momentum L is found to 
have no nonradiative secular change. All spin-spin terms in the secular radiative loss of the energy 
and magnitude of orbital angular momentum are given in terms of L and other constants of the 
motion. Among them, self-interaction spin effects are found, representing the second post-Newtonian 
correction to the 3/2 post-Newtonian order Lense-Thirring approximation. 
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£> ' Neutron-star— black-hole and black-hole— black-hole binaries are among the most promising sources for the earth- 
based gravitational wave observatories, such as the Laser Interferometric Gravitational Wave Observatory LIGO [jlj, 
VIRGO |), GEOjgl and TAMA f§. They are also important sources for the forthcoming Laser Interferometer Space 
I Antenna (LISA) |5||| . It is generally agreed that in a post-Newtonian description of the motion and of the gravitational 
■ radiation of these binary systems a precision of 7/2 post-Newtonian orders has to be reached. At such high orders, due 
to the nonlinear character of the Einstein equations, various multipolar contributions to the gravitational radiation 
7— I . are present. 

Computations based on the general expectation that the orbit of a binary will circularize have reached high post- 
Newtonian orders. However in various physical situations, as described in Refs. 0- |flo|| , the eccentricity of the orbits 
can play a significant role. An equally important feature of the binary system, often neglected, is the effect of the 
individual spins on the motion and on the radiation of the binary. The dynamics under the influence of the spins was 
investigated in Refs. O] - O]. Spin contributions already appear at the 3/2 post-Newtonian order, and exceptionally 
^\ ' at the first post-Newtonian order when angular evolutions are monitored. A complete description should allow for 
Q\ • both spin effects and eccentric motions. 

*"q ' In a recent series of papers Gergely, Perjes, and Vasuth have investigated the influence of the spins in the secular 
evolution and radiation back reaction of binary systems on eccentric orbit up to 3/2 post-Newtonian order. In the 
last paper (ll| of these series, the secular radiative evolution of the energy E and magnitude of the orbital angular 
j^Tj, momentum L were given. Both E and L are constants of the non-radiative motion up to 3/2 post-Newtonian order. 
At this order, among the spin effects, only the spin-orbit interaction is involved in the radiative change of E and L. 
. £h '. The results were in perfect agreement with computations of Rieth and Schafer [jl6f in a different, noncovariant spin 
supplementary condition, with the energy and orbital momentum derived from a different action. 

Also in both the nonradiative and radiative evolution of three angles^] Ki — cos _1 (Si • L), (i = 1,2) and 
7 = cos _1 (Si ■ S2), characterizing the relative orientation of the spin vectors Si and orbital angular momentum L, 
were given. By considering both the spin-orbit and the spin-spin interactions, the nonradiative change of the angles 
receives both first and 3/2 post-Newtonian order contributions, while the radiative changes are of 3/2 post-Newtonian 
order. Remarkably, the radiative change of the spins themselves gave no secular contribution to the radiative evolution 
of the angles. 

Earlier results, pertinent for the physical situations of spinning binaries on eccentric orbit with comparable masses, 
but one of the spins dominating over the second |l7|, a nd a nonspinning particle orbiting about a spinning object (the 
Lense-Thirring approximation ) |]18| (se e also RcfT[|9|), can be obtained as limiting cases of Ref. 15fl. 

The computations in Refs. |15|], |17|, and Jl8| were facilitated by the solution of the radial motion of the binary 
in terms of a generalized true anomaly parameter. Then all relevant integrals presented the remarkable property of 
having the only pole in the origin. Later a systematic analysis pC| has yield the same type of parametrization for a 
wide class of perturbations of the Keplerian motion, including the generic perturbing force of Brumberg [ pl|j22| . 

The main impediment in generalizing the description of Ref. |l5| ] to include all spin-spin effects was that at the 
second post-Newtonian order the magnitude L of the orbital angular momentum, even in the absence of the radiation, 
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An overhat denotes a unit vector. 
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is not a conserved quantity. This feature of the spin-spin type perturbation of the Keplerian motion renders the 
problem even outside the general framework settled in Ref. |B0 |, 

The purpose of this paper is to overcome the above difficulty and to give a solution to the radial motion in terms of 
the generalized true anomaly parametrization, which is valid at the second post-Newtonian order in the presence of 
the spin-spin interaction terms. This task will be completed by a careful analysis of the specific form of the evolution 
of the orbital angular momentum. In Sec. II we study the radial motion. We compute the instantaneous expression 
of L and then we find that it has no secular change. Therefore the angular average value L has a deep significance: 
it has no other change during the whole adiabatic regime of the evolution of the system, then the radiative one. The 
average magnitude L of the orbital angular momentum provides the missing constant of the motion at the second 
post-Newtonian order, and in fact it coincides with the previously employed L up to 3/2 post-Newtonian orders. As 
a by-product, the radial equation of motion is derived. 

Based on the analysis of the radial motion, in Sec. Ill we introduce the generalized true and eccentric anomaly 
parametrizations. The latter is employed in the computation of the radial period, which is found to have the Keplerian 
expression in terms of the energy of the perturbed motion. 

Section IV contains the main results of the paper: the spin-spin contributions to the instantaneous and secular losses 
of the energy and magnitude of orbital angular momentum. The computations are considerably facilitated by the 
advantageous properties of the true anomaly parametrization. Among these losses self-interaction type contributions 
appear. 

In a second post-Newtonian order description of the compact binary there are no other spin terms in the radiative 
losses of energy and magnitude of orbital angular momentum then the spin-spin terms computed in this paper, 
together with the corresponding spin-orbit expressions from Ref. |l5| . We note here that it is equally possible to give 
a description in terms of the time average {L) of the magnitude of orbital angular momentum. However, as will be 
deduced in the Appendix, the relation between (L) and L is cumbersome. The description in terms of L turns out to 
be simpler. We keep the velocity of light c and the gravitational constant G in all expressions. 



II. THE RADIAL MOTION 



The Keplerian motion of a reduced mass particle [i in the gravitational potential generated by the total mass 
to = mi + TO2 is governed by the Lagrangian 



or equivalently, by the acceleration 



/iv 2 Gto/x 

-N = —z 1 , (A-LJ 



Gmr 

sln = o- ■ (2.2) 



Among other dynamical quantities, the Newtonian energy En and magnitude Ln of the Newtonian orbital angular 
momentum are conserved. In terms of the spherical polar coordinates r, 9 and ip they are expressed as 

fiv 2 Gmfi 
En = —z- 



2 r 

^2 +r 2 ( ^ +sin 2„ ,,, «<„(, 
i „• 2 Q -2\ 



-sin 2 (9 (p 2 )] , (2.3) 

r 

L Z N = ^^{0* + sm z e . (2.4) 
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Fr om the first expression of En we find the velocity v. By suitably combining the second expression for En with 
Eq. (2.4) we obtain the equation for the radial motion given by f: 

fi r 
. 2= 2_EN + 2Grn_I^ 
fx r ^ z r z 



When the Keplerian motion is perturbed, Eqs. ( J2.5| ) and (2.6) still hold as identities. However the Newtonian 
expressions En and Ln are not constants of motion any more. The procedure used in Refs. p!q |, pj, and |L8| was to 
express them in terms of the constants of motion E, L, Ki, and 7 of the perturbed motion and of the radius r alone. 
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There is no simple way to fulfill a similar task when the perturbation is due to spin-spin interaction, as will be seen 
in what follows. In fact a milder requirement will be enforced: v 2 and f 2 will be expressed in terms of constants of 
motion and a suitably introduced true anomaly parameter \- 

The motion of the binary under the influence of the spin-spin interactions is governed by the Lagrangian |is)| 



C — £n 
G 



C 



ss 



£ss — 



(Si-S 2 )--(r-S 1 )(r-S 2 ) 



(2.7) 



We have deliberately omitted those first and second post-Newtonian order terms (PN, 2PN) which do not contain 
the spins and the spin-orbit (SO) terms. Their contributions were discussed previously in Refs. |7j] and jig ], and they 
add linearly to the spin-spin perturbative terms computed in this paper. 
An equivalent characterization of the motion is provided by the acceleration 



a = a N + a ss , 



ass 



_3G 

c 2 fj,r 



f {r 2 [r (Si • S a ) + S a (r • S 2 ) + S 2 (r • Si) 



-5r (r ■ S x ) (r ■ S 2 )} 

The perturbed Keplerian motion (2.7) is characterized |l3[ , pi by the conservation of the total energy 

E = E N + Ess 7 Ess = -£.ss 
and of the total angular momentum vector 

J=L+S , 



(2.8) 
(2.9) 

(2.10) 

(2.11) 



where L = Ljv = /i (r x v) is the total orbital angular momentum (there is no spin-spin contribution to the orbital 
angular momentum and the PN, 2PN, SO contributions are not listed here). S = Si + S 2 is the total spin. Both spin 
vectors undergo a precessional motion due to the spin-spin interaction 



Si = 



G 



3("Sj)'-Sj 



x S; 



(2.12) 



where i ^ j. 

From the conservation of the total angular momentum J, Eq. ( 2.11 ) and the precession equations ( 2.12 ) we find 
that the orbital angular momentum evolves as 



3G 



[(r-S 2 )Si + (r-Si)S 2 ] xr 



(2.13) 



The magnitude of the orbital angular momentum changes according to L = L • L, thus 

3G 



L 



-Si 



{ (r-Sa) [MrxSJ 
+ (r-Si)[L-(rxS 2 )]} 



(2.14) 



This change in L is of second post-Newtonian order. Therefore Eq. ( [2.14 ) can be simplified by inserting the previously 
derived Newtonian expressions, Eqs. (2.19) and (2.22) of Ref. |15| : 



r • Si = r sin Ki cos(i/> —i}>i), 
L • (r x Si) = — r sin Ki sin(-0 — ipi) 



(2.15) 
(2.16) 



where the angles ip, ipi are subtended by the node line with the position r a nd th e projections of the spins in the 
plane of the orbit, respectively, see Fig. 1 in Ref. [15[. In order to derive Eq. ( 2.16) the auxiliary relation 



v • Si = f sin Ki cos(V> — ipi ) 



fir 



■ sin Ki sin(-0 — ipi) 



(2.17) 



3 



was employed. The angles tpi are given implicitly in terms of m and 7 by Eq. (2.20) of Ref. [pl| 

S\ sin m cos + S2 sin K2 cos ^2=0 (2.18) 

and by the spherical cosine identity 

cos 7 = coski COSK2 + sinKi sinK2 cos Aip, (2-19) 

where Atp = ip2 — tpi- We introduce the additional notation ip = (tpi + 1^2 )/2. 

In terms of the above angles the variation of the magnitude of orbital angular momentum takes the simple form 

3G 

L = -^r-^Si S2 sin m sin k 2 sin 2(0 — ip) . (2.20) 

c r 

We remark that none of the perturbations not treated here (PN, 2PN, SO) would give contributions to the above 
expression for L. This is to be contrasted with L, which receives also the SO-type contribution [Eq. (2.11) in Ref. 

The change in the magnitude of the orbital angular momentum in the second post-Newtonian order is a feature 
not encountered in the discussion of the PN, 2PN, SO effects. This was the major inconvenience in generalizing the 
description of Ref. |L5| to include the spin-spin contributions. We will see in what follows how to handle this difficulty 
by means of an alternative description based on a properly defined average value L of L. 

The Newtonian limit of a generalized true anomaly parameter x (to be defined later in this paper) has to coincide 
with the Keplerian true anomaly. This feature is enough to introduce it in second post-Newtonian expressions already 
at this stage. As implied by Eq. ( 2.20| ), the magnitude L(x) of the orbital angular momentum at some value of the 
true anomaly \ £ [0, 27r] differs in the second post-Newtonian order from Lq = L (0) : 

f x ■ dt 

L(X) = L + L—d X • (2.21) 
Jo ax 

We compute the integral by inserting Eq. ( 2.20D and the following Newtonian expressions: 

dt ur 2 

T X = -T • ^ 
r = ^ — in V ( 2 - 23 ) 

/i(Gra/i + A cos x) 

^ = X + V'o ■ (2.24) 

The angle t/'o is subtended by the direction of the periastron and the node line, see Fig. 1 in Ref. fjj] and we have 
denoted by L the average value of L(x) over the range (0, 2ir). The quantity A defined by 

i = ^ 2 + ^J 1/2 (2.25) 

is the magnitude of the Laplace- Runge-Lenz vector for a Keplerian motion characterized by E and L. The integration 
yields 

L{x) = L + ^ 2 SxS2 sin k x sin k 2 Q , 

O = 3Gm/x [cos2(?/'o — "0) — cos2(x+V'o — V')] 
+ 4A [0032(^0-^) - cosxcos2(x + V ; o-'0)] 

-2Asinxsin2(x+V'o-'0) • (2.26) 

It is straightforward to verify that L(0) = L(2ir) = L , because all three terms of O independently vanish at x — 0, 27r. 
Therefore there is no change due to the spin-spin interaction in the magnitude of the orbital angular momentum over 
one radial period. 

A second integration over the parameter x from to 2tt gives the angular average of the magnitude of orbital 
angular momentum: 
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where we have introduced the notation 



L=7T L (x)dx = L 



2c 2 Z 3 



(3q = sin Ki sin k 2 cos 2 (-00 — tp) 



(2.27) 



(2.28) 



In terms of L the expression of L(x) reduces to 

Gfi 2 



L(x) = L 



SiS 2 sin^! sinK 2 {2Acos[x+ 2(ip — ip)] 



2c 2 L 3 

+(3GmiJ,+2Acosx) cos 2 (x + 00-0)} 



(2.29) 



A considerab ly sim pler t ask is to find the spin-spin part Ess of the energy. This can be e xpres sed from Eq. (2.1C) 
by use of Eqs. ( 2.15| ) and ( 2.19 ) [the azimuthal angle being eliminated by its expression ( 2.24 )]: 



Ess(r,x) 



G 



■SiS 2 a(x) 



2c 2 r 3 

a(x) = 3 cos ki cos k 2 — cos 7 

— 3 sinKi sinK2 cos 2(x+0q — 0) 



(2.30) 
(2.31) 



Ess is given solely in terms of constants of the motion, radius r an d true anomaly parameter x- Now we are in 
position to write the expressions for v 2 , Eq. (2.5) and for f 2 , Eq. (2.6) in suitable form: 



v 2 = 



2[E-E SS (r,x)} 2Gm 
/i r 

2[E-E SS (r,x)} , 2Gm L( X ) 2 



[l 2 T 2 



(2.32) 
(2.33) 



with L(x) and E ss (r, x) given by Eqs. fl2.29| ) and ( ^30| ). 

As the angles Ki are constant up to the first post-Newtonian order and 7 up to one-half post-Newtonian orde r [see 
Eqs. (2.17 ) in Ref. @], they do not vary in the above expressions. In principle it is possible to rewrite Eqs. ( |2.32| ) 
and (2.33) in terms of constants of motion and r alone. This is because the true anomaly parameter x enters only at 
Newtonian order, in which order the relation x = x( r ) 1S given by Eq. (2.23). We avoid to write such cumbersome 
expressions, as our task will be exactly the converse: to parametrize all expressions in terms of x- 



III. THE TRUE AND ECCENTRIC ANOMALY PARAMETRIZATIONS 



In this section first we introduce the generalized true anomaly parametrization r — r(x) f° r the radial motion 
(2.33). For this purpose we have to find the turning points r max defined by 



r mi „=r(0), r 2 (0) = 0, 



r 



(tt) , f 2 (7r) = . (3.1) 



As from Eq. (2.29) the magnitude of the orbital angular momentum in the turning points is found in the form 
L(0) = L + SL- and L(tt) = L + SL+ with 

SL ± = --^(3Gmn T ^S^o , (3.2) 



the radial equations ( p.l|) evaluated at the turning points have the explicit form: 

n _2E 2Gm _ L(L + 2SL ± ) GS&ap . . 

Here we have denoted a(0) = a(n) = ctQ. We seek for solutions in the form r man . = r± + e±, where e± are small 
corrections to the expressions of the turning points of a Keplerian motion characterized by L and A : 



5 



_ Gra[i ±A_L 2 
T± ~ -2E ~ p{Gmn =F A) ' (3 ' 4) 



The turning points are 

Gmfi ±A Gfi 



-2E 2c 2 AL 2 

p T = au{A^Gmn) + /3 (4A =F 3Gm/x) . (3.5) 



Then we define the generalized true anomaly parameter \ by the relation 

2/1 1 \ / 1 1 



r 



min 1 max / \ 1 min ' max 



cosx , (3.6) 



which yields the expression 

L 2 GfiSiSsA 



fi(Gmfj, + A cos x) 2c 2 AL 2 (Gmfi + A cosx) 2 

with 



(3.7) 



A = A[A 2 (a + 4/3 Q ) + {Gmn) 2 (3a + 10f3o)} 

+ Gm f i[A 2 (3a + llpo) + (Gmfi) 2 (a + 3Po)]cosx ■ (3.8) 

Applications of the true anomaly parametrization ([T?]) will emerge in the next section when the instantaneous 
expressions of radiative losses, represented here by a generic function / will be averaged as 

1 f 27T , dt 



(/) = ^ J f(x)^d X . (3.9) 



In dt/dx — (l/i"){dr/dx) the first factor come s fro m a Taylor series expansion of the radial equation ( 2.33 ) and the 



second is obtained from the derivative of Eq. (3.6): 

1 r 2 sinx- (3.10) 



dr_ _ 1 
~d~X ~ 2 \r. 



m%n ' max 



The parametrization z = exp(i%) has the remarkable property that the integrals (3.£) are given by the residues in the 
origin p0fl . (The only exception under this rule is the time average (L), which is treated in the Appendix.) 

T is the radial period, to be conveniently computed by means of the generalized eccentric anomaly parametrization, 
defined as 



which yields the explicit form 



with 



For the evaluation of the period 



2t {j'rnax ~l~ Tmin) ymax ^min) COS £ , (3.11) 



_ Gray, - Acos£ GpS x S 2 ^ . . 

-2E 2c 2 AL 2 ' [6 ' U) 



S = A(a Q + 4/? ) 4- Gmfx(a + 3/? ) cos£ . (3.13) 



27r dt „ f 2 * Idr 



Jo <K Jo 



the factor 1/r is expressed again from Eq.fl2.33D b y a Taylor series expansion. However this time we have to employ 



also the relations between the Newtonian true and eccentric anomaly parameters 
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Gmucosf — A . J —2EI nL sin£ . 

cos x = -^ a — 1> sm x = ^ r — t ■ ( 3 - 15 ) 

Gm/i — A cos £ Gm/i — A cos £ 



Either the immediate integration or the residue theorem^] yields 



Similarly as in the case of the spin-orbit perturbations, here the period is given by its Keplerian expression (with E 
characterizing the perturbed dynamics). 



IV. SPIN-SPIN CONTRIBUTIONS TO THE INSTANTANEOUS AND SECULAR LOSSES 

The instantaneous radiative losses of the energy E and angular momentum J were given in terms of symmetric 
trace- free (STF) multipole moments by Thornc p3| . To the required order they are 



dE 


G 


dt ~ 


~5c 5 


dJ _ 


2G 


dt ~ 


~5c5 



l(3)jlj(3)jl + 5 j(4\ilk jUlilk + j£j(3)j7j(3)ji 

1 89c 9c 



j(2)j-ij(3)fci + 5 j(3)jimj(4)fcim + 16 j(2).jl j(3)kl 

126c 2 9c 2 



(4.1) 
(4.2) 



where e y fc is the completely antisymmetric Levy-Civita symbol, any number in brackets denotes a corresponding order 
derivative and the expressions of the symmetric trace-free (STF) multipole moments can be derived by use of the 
Blanchet-Damour-Iyer formalism |24| . The Newtonian part of the mass quadrupole moment 1^ and the spin-orbit 
part of the velocity quadrupole moment J J so are given by 



1/ = fi (xix 1 ) 



STF 



J 3l _ 3 ^ 
J SO — o 



Si S2 

mi m,2 



STF 



(4.3) 
(4.4) 



respectively. As estimates of post- Newtonian orders show |14| , these are the only moments needed for evaluation of 
the spin-spin type contribution to radiation losses. 



A. Energy loss 



The Newtonian and spin-spin contribution to the instantaneous radiative loss in the energy is found from Eq. ( |l.l| ): 



dE 
~dt 



G_ 

5c 5 



r(3)ii 
l N 



2G 

5c 5 



In ( a ss)I N > (ajv) + -^J S o ( a w)4o ( a ^) 



(4.5) 



The arguments of time derivatives of the momenta contain the contribution to the acceleration to be inserted in the 
respective terms. The leading order contribution to Eq. (4.5) was given by Pet ers p5| and the S1S2 contribution 
by Kidder |jl4|| , in terms of r,v, S;, and f. Remarkably the last term of Eq. ( |4.5| ) generates also self-interaction 
contributions (containing either or 5|) [p6[. Putting together all these terms, the radiated energy has the expression 



The poles of the integrand (3.14) inside the unit circle of the complex parameter plane w = exp(i£) are in the origin and 

I 1/2 

. The residue at wi vanishes, thus the period (3. If) is given by 



at wx = (emu 2 - y/—2flEL^J / (amy 2 + ^J~^ElPj 
the residue at the origin as T = 2niR(w = 0). 
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dE 
~dt 



8G 3 mV 

15c 5 r 4 
4G 3 m/z 



■(12v 2 - llf 2 ) 



{-171rf[(v • Si) (r • S a ) + (r • Si) (v • S 2 )] 
-r 2 [3(47v 2 - 55f 2 ) (Si ■ S 2 ) - 3 (l68w 2 - 269f 2 ) (r • Si) (r • S 2 ) + 71 (v ■ Si) (v ■ S 2 )]} 

"^S^ E ^ 3r V + 3 ^ + 9f2 ( r ■ S 2 - 6* (r • Si) (v • Si) + r 2 (v • Si ) 2 ] . 



(4.6) 



The summation in the last term is understood over the components of the binary. 

This expression, however is not suitable for averaging by the method described previously. Therefore we insert in 
the Newtonian terms of E q. (|4 .6|) the expressions (2.32) and ( p. 33 ) for v 2 and r 2 , respectively, with L{\) and E(r, x) 
given by Eqs. ( 2.2E) a nd ( |2.30|) . Then we insert in the spi n-spi n terms of Eq. (4.6) the Newtonian expressions for v 2 , 
the Eqs. Q2T5|), (|2.17|) (both with L in place of L), Eq. ( |2.24|) and 



A 



(4.7) 



obtaining the following expression for the instantaneous loss of the energy in terms of the true anomaly parameter \ 
and radius r alone: 



dE_ 
~dt 



dE 

N 



dE\ 

It ) 



N 

8G 3 m 2 



dE\ 
It) 



It) 



SS—self 



S1S2 



. . (2/^EK + 2Gm^r+lLL 2 ) , 
15c°r° 



dE\ 
~dt) 



SS-self 



~ ~ 15c 7 r 8 ^ 



{(6 + sin 2 Ki){8fiEr 2 +8Gm^i 2 r- 3L 2 ) 



+sm 2 K l [4/x^rsinxsin2(x + ipo - 4>i) + (8uEr 2 + 8Gmu 2 r- 5L 2 ) cos2(x + ip - ipi)}} 



dE 
~dt 



2G 3 mS 1 S 2 



SiS 2 



15c 7 (iL 2 t 



a n cos[nx + 2(ipo —ip)] sinKi sinK2 + a± coski COSK2 + as cos 7 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



The coefficients a,k are 



Mu,EL 2 r 2 + 132G 2 m 2 /iV - 52Gm^ 2 L 2 r+ 607L 4 



fll = 4^iAr{33Gmu 2 r- 25L 2 ) , 

a 3 = 4/iAr(llGm/i 2 r+ 25L 2 ) , 

a 4 = L 2 (64:fiEr 2 + 52Gm^ 2 r-465L 2 ) , 

a 5 = L 2 (32/i.Er 2 + 36Gm^ 2 r+ 135L 2 ) . 



(4.12) 



By use of the true anomaly parametrization r(\), Eq. (3.7) we find the energy loss in terms of x alone. When we 
pass to the complex parameter z = exp(i^) the only pole of dE/dt is in the origin. Applying the method described in 
the previous section, the secular change of the energy emerges in terms of constants of motion and the average value 
L of the magnitude of orbital momentum: 



dE 
~dt 

dE\ 
dt I 



dE 



dt ), 



dE\ 
~dt/ 



dE\ 
dt"/. 



N 



dE\ 
~dt) 



SS-self 
dE\ 



SS—self \ I S1S2 

(148£ 2 Z 4 + 732G 2 m 2 u 3 EL + 425G 4 mV) , 

Si 
in, 



G 2 m{-1E^fl 2 
15c 5 L 7 

&m ^lf' 2 E(£) [Ksin'KiCO.Mr., - < , ) + C, Hi + .in 2 ,] 



dt / c 



S-2 



G 2 (-2Eu) 3 / 2 S 1 S 2 
480c 7 L n 



(G3 sinKi sin«2 cos 2(t/>o~ 1 /') + G4 coski cos k 2 + G5 COS7) 



(4.13) 
(4.14) 
(4.15) 
(4.16) 



with the following coefficients Ck 



8 



Ci = -(72E 3 L 6 + 428G 2 m 2 (i 3 E 2 L 4 + AQQG 4 m 4 ^EL 2 + 105G 6 m 6 /i 9 ) , 

G 2 = -4{72E 3 L 6 + 660G 2 m 2 ^E 2 L 4 + 9lOG 4 m 4 n 6 EL 2 + 315G 6 mV) , 

G 3 = -(60744£ 3 Z 6 + 602668G 2 mV 3 ^ 2 i 4 + 819798G 4 m 4 /i 6 ^Z 2 + 266825G 6 mV) , 

G 4 = 4(17064£ 3 Z 6 + 241140G 2 mV 3 -E 2 ^ 4 + 453670G 4 mV£Z 2 + 199899G 6 mV) , 

G 5 = -12(2056£ 3 Z 6 + 28260G 2 m 2 n 3 E 2 L 4 + 52430G 4 m 4 ^EL 2 + 22911G 6 mV) • (4.17) 



B. Loss in the magnitude of orbital angular momentum 



A somewhat more complicated computation yields the similar expressions for the instantaneous and averaged 
radiative losses in L. To complete the task we pass through the following steps. We start from the instantaneous 
radiative loss of J given in Ref . Ej] : 
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We compute the instantaneous radiative loss in L as 

d(3 - S) 
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In Ref. |15||, starting from the Burke-Thorne potential, we have already derived the radiative loss in the spins 
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(4.18) 



(4.19) 



(4.20) 



Qi and being the principal moments of inertia and the angular velocity of the ith spinning axisymmetric body, 
related to the spins by Si — Q^h- These losses appear at the second post-Newtonian order, thus they are not negligible 
in our approach. From Eq. ( 4.2C ) dSi/dt arc found to be perpendicular to the respective spins Si, therefore dSi/dt = 
SidSi/dt. Also in Ref. we have shown that 



dSi 
~~dt 



= 



(4.21) 



This result still holds here, in spite of the fact that the dynamics of the orbital motion and consequently the true 
anomaly parametrization employed for averaging are different. What matters in the derivation of the above result is 
only the Keplerian limit of the parametrizations, which by definiti on, co incide. 

Therefore in the secular loss (dL/dt) only the first term of Eq. ( 4.19 ) is relevant, which will be computed in what 
follows. As PN , SO, and 2PN contributions are not considered in our t reatme nt, w e insert Ljv= L in the Newtonian 
term of ( 4.18 ), together with the previo usly derived expressions ( 2.32 ) and ( 2.33 ) for v 2 and f 2 , respectively. The 



spin-spin perturbation term of Eq. (4.18) can be easily dealt with by use of the Cartesian coordinates (x — r cost/', y — 
r sin?/;, z = 0). Then to the r equire d, le ading order L = (0, 0, 1) (see also Fig.l. of Ref. (l5)). We employ in addition 
the Keplerian relations (2.22) and (2.24 ) which imply 



(4.22) 



A lengthy but straightforward computation yields the following implicit function of the parameter x f° r L-d3 /dt: 

(4.23) 
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4G 2 S 1 $ f x 3 ■> - ] 

= T1~UT1 \ lli ^ n C0S [ n X + 2(V>o -"0)] siiiKi sinK 2 + 64 coski cosk 2 + 65 cos 7 > , (4.26) 
SiSa oc fi L r y n=1 J 

with the coefficients given by 

61 = 3fiAr(-2Gmfi 3 Er 3 + f iEL 2 r 2 + 9Gm[i 2 L 2 r - 8L 4 ) , 

b 2 = 3(-2G 2 m 2 n 5 Er 4 - 22fiEL 4 r 2 + 9G 2 m 2 fi 4 L 2 r 2 - 10Gm/x 2 £ 4 r + 19L 6 ) , 

63 = fiAr{-2Gmfi 3 Er 3 - 3^EL 2 r 2 + 9Gm^i 2 L 2 r + 2AL 4 ) , 

64 = Z 4 (54^£V 2 + 23Gmfi?r - 45L 2 ) , 

65 = -3Z 4 (6^£V 2 + 2Gm[i 2 r - 5L 2 ) . (4.27) 



We would like to stress the appea rance of the self-interaction terms ( 4.25 ) in the instantaneous loss of L either. They 
originate in the last term of Eq. ( 4.1S| ) . 

A similar averaging procedure as in the case of the energy loss, gives here 



dL\ _ /dL\ I dL\ /dL\ 
dt J \dt/ N \dt/ ss _ self \dt/ SxS2 ' 

dL\ 4G 2 m(-2fiE) 3 / 2 M - a , - 2 3 . 

-dt) N = h£^ {UEL + 15G mV) ' 

dL\ G 2 mii{-2iiE) 3 ' 2 „ v-^ / Si \ 2 „ . 2 N 

<2L\ G 2 (-2fiE) 3 / 2 S 1 S 2rT ^ . . tn „ - 

— ) = =r=- Lf 2 sinKi sin(t 2 cos2(^ — V) + ^3 cos K i cos/c 2 + cos 7 

I SiS 2 lOc'i 8 

where the coefficients are given below 



(4.28) 
(4.29) 
(4.30) 
(4.31) 



D x = \2E 2 L 4 + 60G 2 m 2 n 3 EL 2 + 35G 4 mV , 

D 2 = 12(62£ 2 Z 4 + 2\\G 2 m 2 p?EL 2 + 90G 4 m 4 ^ 6 ) , 

D 3 = -(1524£ 2 Z 4 + 7260G 2 m 2 n 3 EL 2 + 4865G 4 mV) , 

L> 4 = 24(22£ 2 Z 4 + 105G 2 m 2 n 3 EL 2 + 70G 4 m 4 ^ 6 ) . (4.32) 



We emphasize that self- interaction terms (4.15) and (4.30) are present in both the secular losses of E and L. 



V. CONCLUDING REMARKS 



The main result of this paper consists in the derivation of the the second post-Newtonian order spin-spin contribution 
to the radiation back reaction on the energy and magnitude of orbital angular momentum of the coalescing binary 
system. In order to complete the task, a solution r(x) of the quasi-Keplerian radial motion of the binary, with the 
spin-spin interacti on as t he p erturbation, was deriv e d in terms of the true anomaly parametrization. The secular 
losses of E, Eqs. ( 4.13| )-( |lTf ) and of L, Eqs. ( 4.28 )- (|4.31 ), in complement with the corresponding spin-orbit terms 
from Eqs. (4.1), (4.2) of Ref. jl5| (with the replacements L — > Z and A — * A, which are allowed up to the 3/2 post- 
Newtonian order), represent the total spin- induced back reaction of the radiation up to the second post-Newtonian 
order. Among the loss terms due to the spin-spin interaction, terms originating in the self-interaction of the spins were 
found. While all other spin-spin terms vanish in the one-spin limit S 2 = 0, (half of) these terms survive, representing 
the second post- Newtonian correction to the losses derived earlier in the Lense-Thirring approximation [ p"8"| . 

In the nonradiative description of the angles Ki and 7 in principle not all second post-Newtonian order spin terms 
can be found from the Lagrangian (2.7). Similarly as the leading spin-orbit and spin-spin terms in the spin precession 
equations (exceptionally) are of first and 3/2 post-Newtonian orders, respectively jL5|| , any higher order description of 
the evolution of the spinning binary then Eq. ( |2.7| ) could give new contributions even at the second post- Newtonian 
order. Such a description is in the process of completion pq. On the other hand, the situation is more sophisticated 
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when it comes to the radiative evolution of the angles Ki and 7. New terms beyond those already given in Ref. M] can 
arise in two ways. The first type of these terms would come from corrections to the Burke-Thorne potential [^7| . It was 
already shown 1 29 1 that the first correction is one post-Newtonian order higher than the potential itself, therefore any 
correction from the radiative loss of the spins would contribute only at 5/2 post-Newtonian order (above the leading 
radiation term in the orbital angular momentum loss). However new terms in the radiative angular evolution appear 
at the second post- Newtonian order due to (dJ /dt)ss- A computation in progress reveals a complicated structure of 
these terms. We defer the study of the angular evolutions at the second post-Newtonian order to a subsequent work. 
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APPENDIX A: THE RELATION BETWEEN TIME AND ANGULAR AVERAGES OF L 

We start from the expression ( f2~29| ) for L(x). The time average is obtained as 



Here T is given by Eq. (3.16) and dt/dx is derived from the radial equa tion as explained in Sec. III. When expressed 
in terms of the complex parameter z = exp(ix), the integrand in Eq. (Al) has two poles inside the unit circle: one 
in the origin zq — and the other at 



Zl 



-Wl 



Gray 2 - ^ -2^EL 2 
Gmp 2 + yJ-I^EL 2 



1/2 



(A2) 



From the residue theorem, {L) is given as the sum of the residues at zq and z\ multiplied with 2iri/T. 
Thus the desired relation between the angular average L and the time average (L) emerges as 



(L)=L 



2c 2 A 2 L 3 E 2 



5*i 5*2 sinKi sin k 2 cos2(?/;o — ip) 



(A3) 



with the coefficients E\^. given by 



E x = 2(-2 ijlE) 1 ' ' 2 L[A 6 - 15G 2 7n 2 ^ 2 A i + 32GSn i ^A 2 - 16G 6 mV] 

+Gmn 2 [~llA 6 + 58G 2 m 2 n 2 A 4 - 80G 4 m 4 n 4 A 2 + 32G 6 mV] 
E 2 = 4Gm{-2nE) 1/2 L[A 2 - 2G 2 m 2 ^ 2 ] + A A - 8G 2 m 2 fi 2 A 2 + 8G 4 m 4 /i 4 



(A4) 



Whenever needed, a ll ex pression s can be rewritten in terms of (L) by use of Eq. (A3). However, by inserting L 
expressed fro m Eq . ( |A3| ) in Eq. Q2.29D we obtain a considerably longer expression for L(\) m terms of (L) than the 
original Eq. ( 2.29 ). This is an indication that the description in terms of the angular average L is better suited for 
the problem. 
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